Hadron production in deuteron-gold collisions and nuclear parton distributions 
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We calculate nuclear modification factors RdAu, central-to-peripheral ratios, Rcp, and pseudo- 
rapidity asymmetries YAsym in deuteron-gold collisions at y/s = 200 GeV in the framework of 
leading-order (LO) perturbative Quantum Chromodynamics. We use the Eskola-Kolhinen-Salgado 
(EKS), the Frankfurt-Guzey-Strikman (FGS) and the Hirai-Kumano-Nagai (HKN) nuclear parton 
distribution functions and the Albino-Kramer-Kniehl (AKK) fragmentation functions in our calcu- 
lations. Results are compared to experimental data from the BRAHMS and STAR collaborations. 

PACS numbers; 24.85.-|-p,25.30.Dh,25.75.-q 



I. INTRODUCTION 

In the rapidly-developing field of relativistic nuclear 
collision physics, questions related to the distribution 
of partons in nucleons and nuclei are of great current 
interest. The modification of the parton distribution 
functions (PDFs) in the nuclear environment has at- 
tracted growing attention ever since the pioneering EMC 
experiment [U, At the Relativistic Heavy Ion Col- 
lider (RHIC), we are well-positioned to study the nu- 
clear PDFs (nPDFs) in a center-of-mass energy and 
transverse-momentum domain where perturbative Quan- 
tum Chromodynamics (pQCD) is expected to work well. 
Thus, hadron production in nuclear collision experiments 
at RHIC should provide information on how such non- 
perturbative ingredients of pQCD as the parton distribu- 
tion functions are modified by the presence of the nuclear 
medium. 

However, the description of nuclear collisions based on 
pQCD is a complicated task. Much of the complication 
derives from the intrinsically complex nature of the nu- 
clear environment in collisions of heavy nuclei. In ad- 
dition to the initial-state modification of the PDFs, the 
cross section of high-p^ hadron production is infiuenced 
by final-state effects (such as jet energy loss) and a com- 
plicated geometry. To better understand the physics of 
pQCD in the nuclear environment, it is highly desirable 
to disentangle the different nuclear phenomena affecting 
high-pT [pt ^ 2 GeV/c) hadron production. These phe- 
nomena (not present in proton-proton collisions) mani- 
fest themselves in e.g. the measured nuclear modification 
factors. 

Deuteron-gold (d-|-Au) collisions provide a good com- 
promise and testing ground for these purposes. The 
deuteron, being the simplest "real" nucleus, affords a 
complexity higher than a proton, but much less than that 
of a typical heavy nucleus like gold. Thus, a good under- 
standing of d-f-Au collisions is invaluable in elucidating 
the added complexity associated with collisions of more 
complex heavy nuclei. Accordingly, d-|-Au collisions have 
been extensively studied at RHIC (see e.g. 0, 13, H, 111 ) ■ 
A new feature, offered by nonidentical colliding beams 
like d+Au, is the pseudorapidity asymmetry, examined 



in some detail recently by the STAR collaboration^. 

Any pQCD calculation involves, in addition to par- 
tonic differential cross sections, parton distribution func- 
tions (PDFs) and fragmentation functions (FFs) to con- 
nect to the observable level. The latter ingredients are 
non-perturbative, but universal in the absence of the nu- 
clear environment. Typical pQCD-based calculations of 
nucleus-nucleus collisions use modified (nuclear) PDFs 
and deal with issues like jet energy loss (jet quenching) 
and collision geometry. In this work, we focus on the 
PDFs (and of course can not avoid treating the geome- 
try of the collision). Even in the proton (nucleon), our 
knowledge of the PDFs is naturally limited; the nPDFs 
arc much less well known in the wide range of momentum 
fraction x needed for reliable calculations. The nuclear 
gluon distribution, in particular, is poorly constrained. 
The uncertainties in the nPDFs directly affect the accu- 
racy of pQCD calculations. It is therefore important not 
to rely on a single nPDF parametrization. Calculations 
utilizing different parameterizations offer a useful check 
on the performances of the different nPDFs in describing 
relevant observables. 

In the present work we compute nuclear modifications 
expressed in terms of the ratios RdAu (nuclear modifica- 
tion factor, see eq. ([8])) and Rcp (central-to-peripheral 
ratio, see eq. pT|) ). We also calculate the pseudorapid- 
ity asymmetry YAsym (see eq. (|12|)) in certain pseudora- 
pidity intervals. We focus attention on the phenomenon 
of nuclear shadowing, the difference between PDFs and 
nPDFs, leaving aside possible additional effects like jet 
quenching or intrinsic parton transverse momentum and 
its broadening in nuclear collisions. In this way we estab- 
lish a minimalist base line for experimental comparisons. 

An earlier study along these lines [1] used the 
Eskola-Kolhinen-Salgado (EKS)[1] and Frankfurt-Guzey- 
Strikman (FGS) [IQ] parameterizations of nuclear shad- 
owing. In Refs. [9|, [lO| the basic object is the shadow- 
ing function, which encodes the relevant nuclear infor- 
mation. The nPDFs are then a product of the shad- 
owing function and nucleon PDFs. In our present study 
we are particularly interested in the Hirai-Kumano-Nagai 
(HKN) parameterization 11], which also gives the par- 
ton distribution functions for the deuteron, unlike other 
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parameterizations in which the deuteron is not shad- 
owed. (Deuteron shadowing is of course a small effect 
compared to shadowing in heavy nuclei.) A short re- 
view of the available nuclear parton distribution func- 
tions and their differences can be found in Ref. An- 
other point of departure in the present study from Ref. [1| 
is the use here of the Albino-Kramer-Kniehl (AKK) frag- 
mentation functions [Tsj . These updated fragmenta- 
tion functions are only available in next-to-leading or- 
der (NLO) and incorporate new experimental informa- 
tion from the OPAL Collaboration, including light quark 
tagging probabilities 14]. Thus AKK is expected to of- 
fer a better description of the fragmentation process. For 
calculations involving the EKS and FGS shadowing func- 
tions we need the nucleon parton distributions. We use 
the MRST2001 leading order (LO) PDFs[ill for consis- 
tency with HKN, where the underlying nucleon PDFs are 
the MRST2001 LO PDFs. 

The paper is organized as follows: in Sec. |TT]we review 
the basic formalism of LO pQCD as applied to d-|-Au 
collisions. This Section also includes the definitions of 
the nuclear modification factor, the central-to-peripheral 
ratio, and the pscudorapidity asymmetry, and a brief 
review of the available experimental data. We present 
the results of our calculation in Sec. IIIIl and conclude in 
SecHVl 



II. BASIC FORMALISM 

The invariant cross section for the d-|-Au h+X re- 
action, with respect to pscudorapidity 77 and transverse 
momentum px can be written as 



Ma 

drj(PpT 



^ = ^ Jd^b d^s dxadxbdzc tdis) tAu{\b - s\) 

abed 
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-s5{s + t + u) ,(1) 



where Xa and Xb are parton momentum fractions in 
deuteron and gold, respectively, and Zc is the frac- 
tion of the parton momentum carried by the final-state 
hadron h. The factorization and fragmentation scales are 
Q and Qf, respectively. Here, 



tA{s) = dzpA{s,z) 



(2) 



is the Glauber thickness function of nucleus A, with the 
nuclear density distribution, pa{s, z) subject to the nor- 
malization condition 



J d^s dzpA{s, z) = A 



(3) 



The quantity d(j(ab cd) /di in eq. |T]) represents 
the perturbatively calculable partonic cross section, and 



Dh/cizcQf'^) stands for the fragmentation function of 
parton c to produce hadron h, evaluated at momentum 
fraction Zc and fragmentation scale Qf. Using the i5- 
function in eq. ([T]), the integration over Zc can be car- 
ried out explicitly. Integration limits over Xa and Xb are 
then (xamm, 1) and (xbrmnixa), 1) respectively. Note that 
function of Xa- In addition, Zc is also a function 
of both Xa and Xb- 

In the present study, we are primarily concerned with 
Fa/Ai^i z), the nuclear parton distribution function 
(nPDF) for nucleus A. In light of the nuclear mod- 
ifications discussed in Sec. |T1 it is natural to assume 
that the nPDF depends on the location of the parton 
in the nucleus, (s, z) (or at least on its position rel- 
ative to the beam axis s). To connect this "inhomo- 
geneous" nPDF to the geometry-independent (homoge- 
neous) nPDF ^/yi(x, Q^), the normalization condition 



d^s dz pa{s, z)Fa/A{x,Q^ ,s, z) = Ta/A{x,Q'^) (4) 



should be satisfied. 

In the EKS (3, FGS[l3, and HIJING[i11 parameteri- 
zations the (homogeneous) shadowing function S{x, Q^) 
is introduced, and the nPDF is written as 



J^a/AiX; Q^) S{x, Q^) fa/Nix, Q^) 



(5) 



where fa/N{x,Q'^) is the PDF of the nucleon, which can 
be expressed as 

fa/N{x,Q^) - jfa/pix^Q^) + (1 - j).fa/nix,Q^) , (6) 

with fa/p{x, Q^) [fa/n{x, Q^)] being the proton [neutron] 
parton distribution function as a function of Bjorken x 
and factorization scale Q. The HKN parameterization is 
already given in terms of nPDFs. 

Especially at small values of x, where coherence ef- 
fects are important, scaling with the thickness function 
appears to be more physical than with the local density. 
Therefore, in this work we assume that shadowing is pro- 
portional to the thickness function ^ (i.e. we adopt the 
second option discussed in Ref. [8|). 

The connection between the inhomogeneous and 
homogeneous shadowing functions, Sa{x,Q'^,s,z) and 
S{x,Q^), respectively, can be written as 



SAix,Q^,s,z) ^l + N[S{x,Q^) 



JdzpA{s,z) 
^JdzpA{0,z) ' 



with N a normalization constant. In other words, the 
deviation of the inhomogeneous shadowing function from 
unity is proportional to the deviation from unity of the 
homogeneous shadowing function and proportional to the 
thickness function. While inhomogeneous FGS nPDFs 
are available, for consistency and ease of comparison with 
EKS and HKN, we apply the homogeneous FGS nPDFs. 

We obtain the density distribution of the deuteron 
from the Hulthen wave function[17] (as in Ref. 
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while a Woods-Saxon density distribution is used for 
gold with parameters from Ref. [l^. Since the Nijmegen 
deuteron wave function[2^, which we also applied, gives 
similar results to the Hulthen wave function, we report 
only the calculations using the Hulthen wave function 
here. 

We fix the scales as Q — Qf — pr, where px is the 
final hadronic transverse momentum. We also carried 
out calculations with the scales Q = px/zc, Qf — pr- 
Results with the latter choice do not differ significantly 
from those obtained by having both scales fixed at px- 
The partonic differential cross sections, da{ab cd)/dt 
were evaluated at leading order (LO). We note that, if a 
K factor was used to approximate the effects of higher 
orders, these effects would cancel in the ratios calculated 
in the present study. For the fragmentation functions we 
use the AKK set[l^ throughout. 



Nuclear Modification Factors 



The d+Au nuclear modification factor, RdAu is defined 



The nuclear modification factors RdAu have been mea- 
sured at several pseudorapidities by the BRAHMS Col- 
laboration, and are presented at I77I < 0.2(77 = 0), 
0.8 < ?/ < 1.2(77 = < V < 2.35 (?7 = 2.2), and 

2.9 < r? < 3.5 {-q = 3.2)[||. At small rapidities {-q = 0, 1), 
the data are given for the sum of charged hadrons, while 
at forward rapidities negatively charged hadron data are 
available. The AKK fragmentation functions are for av- 
erage charged hadrons. It should be remembered that the 
calculated results for all rapidities are therefore for the 
average of charged hadrons. Furthermore, the BRAHMS 
data are given in three centrality classes: central (0- 
20)%, semicentral (30-50)%, and peripheral (60-80)% 
(as a percentile of the geometric cross section). With 
a Woods-Saxon density for gold and the Hulthen wave 
function for the deuteron, a Glauber calculation of TdAu 
relates these classes to impact parameter intervals asQ 
< < 3.81 fm for central, 4.66 < b < 6.01 fm for semi- 
central, and 6.59 < 6 < 7.74 fm for peripheral. We will 
use the notation Rcp to refer to the central-to-peripheral 
ratio, while Rsp will be used to denote the semicentral- 
to-peripheral ratio in the following. 



as 



RdAuiPr) 



1 'i'^dAu 



da, 



pp 



(Nbtn) drj d pt / '^V (Ppi 



(8) 



where the average number of binary collisions, (Nbin) in 
the various impact-parameter bins is given by 



(Nbin) = {(^NN TdAu{b)) 



(9) 



Here ct]^^ is the inelastic nucleon-nucleon cross section, 
and 

TdAu{b) = Jdh td{s) tAu{\b- ^) (10) 

represents the deuteron-gold nuclear overlap function. 
The nuclear modification factor RdAu is thus just the 
ratio of the d-f Au and proton-proton (pp) cross sec- 
tions, normalized by the average number of binary colli- 
sions, {Nbin)- 



B. Central-to-Periplieral Ratios 

A related ratio, which dispenses with the need for a 
reference pp cross section and uses information from the 
same experiment in numerator and denominator, thus 
canceling most systematic errors, is the central to pe- 
ripheral ratio defined as 



Rcp{pi 



1 



dcT^^ 

^'^dAu 



1 



drrf^P 
'^'^dAu 



{Nbin)c drj d^pT I {Nbi7i)p drj d' 



Pt 



(11) 

where (Nbin) is as defined above. The label C stands 
for the central event class, while P denotes the periph- 
eral class. The centrality classes are chosen according to 
centrality cuts on the experimental data. 



C. Pseudorapidity Asymmetry 

As the mechanisms for hadron production in d+Au 
collisions may be different at forward rapidities (deuteron 
side) and backward rapidities (gold side), it is of inter- 
est to study ratios of particle yields between a given 
rapidity value and its negative in these collisions. The 
STAR Collaboration has recently measured pseudorapid- 
ity asymmetries Q, defined as 



YAsym = ^"b'' (Au-Side) / ^"^^^'^ (d-sidc) 



sym — , ,2 

dr] d PT 



dr] d PT 



(12) 



in d+Au collisions for several identified hadron species 
and total charged hadrons in the pseudorapidity inter- 
vals \ri\ < 0.5 and 0.5 < |?7| < 1.0. Rapidity asymmetries 
with the backward/forward ratio above unity for trans- 
verse momenta up to « 5 GeV/c are observed for charged 
pion, proton-|-anti-proton, and total charged hadron pro- 
duction in both rapidity regions. We want to see if dif- 
ferent nPDFs give significantly different rapidity asym- 
metries for the various hadron species. 



III. RESULTS 

Before presenting the results from the three different 
nPDFs, we note that while EKS and HKN are similar in 
the sense that they are global fits to experimental data, 
FGS is relying on Gribov theory for diffractive deep in- 
elastic scattering to derive nPDFs. Gluon shadowing is 
much stronger in FGS than in EKS and HKN. Because 
nuclear gluon distributions arc poorly constrained exper- 
imentally, there are significant differences between the 
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EKS and HKN gluon shadowing. On the other hand, 
since the Gribov formahsm is not capable to predict va- 
lence quark shadowing, FGS uses EKS shadowing for 
valence quarks, and all three parameterizations are in 
agreement for valence quark shadowing. For sea quarks, 
EKS and HKN are in good agreement at 0.01 ^ x < 0.1. 
At a; > 0.2 and in the small-x region x < 0.01 there 
are substantial deviations. The FGS predicts more nu- 
clear shadowing for sea quarks than EKS and HKN. Due 
to the x-integration, different nuclear effects [shadowing 
(x < 0.1, depletion), antishadowing (0.1 ^ x < 0.3, en- 
hancement), EMC effect (0.3 < a; < 0.7, depletion), etc] 
are superimposed, thus it is difficult to isolate these ef- 
fects. 

Since pQGD calculations are generally not reliable at 
low pt, we do not wish to push our calculations below 
Pt = 1-5 GeV/c. With our scale choice, this corresponds 
to a minimum (Qj) of 2.25 GeV^, while EKS, FGS, 
and HKN give 2.25,4.0, and 1.0 GeV^ for minimum Q^, 
respectively. The minimum for the AKK fragmenta- 
tion functions is given as 2.0 GeV^. 

A. Nuclear Modification Factors 

We have calculated the minimum bias RdAu for total 
charged hadron production at the BRAHMS pseudora- 
pidities with the three nuclear parton distribution func- 
tions considered in this study. In the case of the FGS 
nPDFs we use the strong gluon shadowing for the gold 
nucleus. We have only calculated the px distributions 
for final-state total charged hadrons in the present study, 
since we are mainly interested in the performance of the 
different nPDFs. More detail for different hadron species 
and fractional contributions from quarks, antiquarks, and 
gluons for EKS and FGS nPDFS can be found in [|. 
Here we limit our discussion of nuclear modification fac- 
tors to the effect of the different nPDFs. We make two 
general observations: 

i) At midrapidity (small 77), processes initiated by 
both gluons and quarks are important. At forward 
rapidities (large 77), xt, the parton momentum frac- 
tion in gold becomes small and gluon-initiated pro- 
cesses become dominant. 

ii) In both, the data and the calculations, RdAu de- 
creases systematically with increasing 77. This re- 
flects the increasing role of shadowing since smaller 
values of Xb are probed at forward rapidities. 

The results of our calculations are displayed in Fig. [l] 
together with the experimental data. At 77 = and 
Pt = 1-5 GeV/c, typical Xmin values are significantly be- 
low 0.1. Thus there is substantial contribution from the 
shadowing region, and all three nPDFs predict RdAu < 1, 
with the HKN being the lowest. Due to their steeper rise 
of RdAu with Pt, the EKS and FGS nPDFs appear to de- 
scribe the data better at around pT — 3 GeV/c than the 



HKN nPDFs. At pt = 4 GeV/c, Xbrmn is above 0.5 for 
some values of Xa- There is thus more contribution from 
the antishadowing region, and consequently RdAu > 1, 
with the HKN stiU being the lowest. At pT = 8 GeV/c, 
we have major contributions from both antishadowing 
and the EMC effect. Thus, while RdAu is still > 1, the 
trend is towards 1 for both EKS and FGS. The HKN 
parameterization predicts a higher value at around 1.1. 

The behavior is similar at = 1. At pT = 1-5 GeV/c 
both Xmin values are small, and thus RdAu < 1- Around 
this Pt the HKN is identical with the FGS. At pT = 
4 GeV/c, Xbmin ranges higher, with RdAu > 1- The HKN 
is slightly below both EKS and FGS in this region. At 
Pt = 8 GeV/c, RdAu > 1 with substantial contributions 
from both antishadowing and the EMC effect. The HKN 
is practically the same as EKS and slightly higher than 
FGS. All three are in good agreement with data except 
at low Pt, where the data are more suppressed than the 
calculated results. 

The effect of increasing ry is already apparent at 77 = 
2.2, where both FGS and HKN are below unity for all 
Pt considered. The EKS is practically 1 at mid-p^ and 
falls below 1 at high pT- The major contribution is from 
shadowing with the resultant RdAu < 1- Note that at 
low Pt the HKN is similar to the EKS, while at higher 
Pt it is lower than both EKS and FGS. The agreement 
with data is reasonable. 

At r; = 3.2 the dominant contribution is again from 
shadowing with RdAu < 1 for all pT- At low pT the FGS 
describes the data best while at high pT it is similar to the 
EKS. This may be due to the stronger gluon shadowing in 
the FGS. The HKN is less suppressed at low pT than both 
EKS and FGS, while at higher pT it is more suppressed 
than the others.. 



B. Central-to-Peripheral Ratios 

The results of our central-to-peripheral ratio (Rcp) 
calculation are displayed in Fig. [2l together with the 
BRAHMS data. At 77 = 0, where the data indicate an 
Rcp > 1 at 2 < Pt ^ 4 GeV/c, the calculated results 
are below the data, while at both rj = 2.2 and 77 = 3.2 
the data show significant suppression, with the calcula- 
tions giving Rcp close to unity. At 7; = 1 the calcula- 
tion can be said to be in reasonable agreement with the 
data. Although the calculation exhibits the trend toward 
increasing suppression in the data as rj increases, the cal- 
culated variation with 77 is much smaller than the one 
shown by the data. This shortcoming becomes increas- 
ingly evident at forward rapidities. Fig. [3] shows the cal- 
culated semicentral-to-peripheral ratio, Rsp, using the 
three nPDFs, together with the BRAHMS data. The 
situation here mirrors that of the central-to-peripheral 
ratios. The results underpredict the data at rj — 0, but 
overpredict at forward rapidities. The degree of suppres- 
sion at forward rapidities is not as severe as in the central- 
to-peripheral ratios. In fact, due to the large error bars, 
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FIG. 1: (Color Online) Minimum bias nuclear modification factor RdAu for total charged hadron production at different 
pseudorapidities, using EKS (solid), FGS (dashed), and HKN (dot-dashed) nPDFs. Solid triangles denote BRAHMS data. In 
the top panels the data are for average charged hadrons, while for the bottom panels only negative hadrons are measured. The 
solid error bars are systematic errors, while dashed error bars represent statistical errors. For clarity, only statistical errors are 
displayed in the two upper panels, where the systematic errors are rather large. 



there is reasonable agreement with data at pT around 
4 GeV/c. 

Both Rcp and Rsp are geometry-dependent, thus the 
assumed spatial dependence of shadowing is important. 
We have checked that shadowing proportional to the lo- 
cal density gives worse agreement with data than shad- 
owing proportional to the thickness function. A varia- 
tion of the thickness function dependence where we used 



higher powers gives better agreement at 77 = but still 
overpredicts significantly at other rapidities. With these 
shadowing parameterizations, a more radical spatial de- 
pendence is needed to describe the data for both ratios. 
Another factor that may be responsible is a too-weak 
x-dcpendence of the available shadowing parameteriza- 
tions. 
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FIG. 2: (Color Online) Central-to-peripheral ratio Rcp for total charged hadron production at different pseudorapidities using 
EKS (solid), FGS (dashed), and HKN (dot-dashed) nPDFs. Solid triangles denote BRAHMS data. The solid error bars are 
systematic errors, while dashed error bars represent statistical errors. 



C. Pseudorapidity Asymmetry 

In d+Au collisions at small rapidities, particle pro- 
duction may include contributions from gold-side par- 
tons that may have been modified by nuclear effects and 
from deuteron-side partons that have experienced multi- 
ple scatterings while traversing the gold nucleus [7|. It 
should be kept in mind that the latter effect is not in- 
cluded in the present calculations. 

We have calculated Yasym for charged pions (7r+ + tt^ ) , 



charged kaons, protons-|-antiprotons {p + p), and total 
charged hadrons + h~). Below 2 GeV/c, data are 
available for total charged hadrons only. Above 2 GeV/c, 
separated data exist for 7r+ -I- tt" and for p + p. A ben- 
efit of using a ratio is that the systematic errors largely 
cancel and are < 5% for both pions and protons, < 3% 
for charged hadrons 7]. For pT S 4-2 GeV/c the errors 
are dominantly systematic and thus the statistical errors 
are not displayed. At higher pr statistical errors tend 
to become dominant. Since the asymmetry is a ratio of 
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FIG. 3: (Color Online) Semicentral-to-peripheral ratio Rsp for total charged hadron production at different pseudorapidities 
using EKS (solid), FGS (dashed), and HKN (dot-dashed) nPDFs. Solid triangles denote BRAHMS data. The solid error bars 
are systematic errors, while dashed error bars represent statistical errors. 



the yields in two different rapidity intervals, the respec- 
tive Xfc-distributions are mostly responsible for the global 
trends observed in the calculations. 

Fig. m shows the result of our pseudorapidity asymme- 
try ratio calculations for the interval I77I < 0.5, together 
with the STAR dataQ- We plot the calculated results 
with the EKS (sohd line), FGS (dashed), and HKN (dot- 
dashed) nuclear parton distributions. Around 77 = one 
expects the degree of asymmetry to be small. This is 
borne out by both data and calculation. At very low 



Pt, there is more contribution from the shadowing region 
for the positive rapidity (deuteron-side) than the nega- 
tive rapidity (gold-side). Thus the ratio is expected to 
be above unity. For px > 3 GeV/c, the xi, distributions 
are similar, with less contribution from the shadowing re- 
gion as Pt increases. Thus the asymmetry is not far from 
unity. The EKS and FGS nPDFs give similar results for 
all hadronic species. The HKN nPDFs yield a different 
curvature and pseudorapidity asymmetries that remain 
above unity for a wider range of transverse momenta. 
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For the total charged hadrons ratio, the calculation is 
below the data for pt ^ 5 GeV/c. At the present level 
of combined experimental and theoretical uncertainties, 
we can claim that Yasym is reasonably reproduced for 
TT"^ + TT^ and p + p, in particular at large transverse mo- 
menta. The EKS and FGS predictions for charged kaons 
have a small maximum at px ~ 2.2 GeV/c. 

Figure [5] shows the result of our pseudorapidity asym- 
metry ratio calculations for the interval 0.5 < \ri\ < 1.0. 
At low pt, the situation is similar to {rjl < 0.5. Thus the 
asymmetry is above unity. However, at > 3 GeV/c, 
the Xb distributions start to become significantly differ- 
ent. While there are still some contributions from the 
shadowing region for the positive rapidity even up to the 
highest Pt, the negative rapidity yield has no shadow- 
ing contribution for Pt ^ S> GeV/c. Thus one expects 
the asymmetry to be more substantial than for |?7| < 0.5. 
This is borne out by the calculations. As in |?7| < 0.5, 
EKS and FGS give very similar results. The HKN nPDFs 
also behave similarly to what was seen at |?7| < 0.5, the 
difference between EKS and FGS on the one side, and 
HKN on the other, becoming more pronounced. All cal- 
culations underpredict the data for h^ + h~ . For tt+H-tt^ 
Build p + p, the calculation agrees with the data within er- 
ror for large transverse momenta. It will be interesting 
to see the data for charged kaons, when they become 
available [2l|. 

IV. CONCLUSION 

We have calculated nuclear modification factors RdAu 
in deuteron-gold collisions at -/s — 200 GeV using three 
nuclear shadowing parameterizations: the EKS, FGS, 
and HKN nuclear parton distribution functions. We have 
also calculated a related modification factor, the central- 



to-peripheral ratio, Rcp, using the impact parameter 
representation of the BRAHMS centrality classes. We 
have used the AKK fragmentation functions throughout 
in our calculations. 

All three nPDFs give similar results for the minimum 
bias RdAu at forward rapidities. For not too large rapidi- 
ties and relatively low pt, the HKN values are smaller 
than those of the EKS and FGS. The calculated results 
for all three nPDFs approximately agree with the data. 

The Rcp results do not describe the data well over- 
all. At \ri\ < 0.2 the calculated results underestimate the 
data, while at the most forward rapidities the data are 
much below the calculated results. Only at r/ = 1 can 
we claim approximate agreement with data, but this ap- 
pears to be coincidental, judging from the 77-dependence 
of the data and the calculated results. This situation 
may be due to poor knowledge of the spatial dependence 
of nuclear shadowing. It may also signify a rapidity de- 
pendence of shadowing. 

For all hadronic species the calculated pseudorapidity 
asymmetry underestimates the data for pt ^ 5 GeV/c. 
At higher pT, EKS and FGS results tend to fall lower 
than unity, while HKN yields values above one. Present 
uncertainties in the experimental data do not allow a 
distinction between these scenarios. We expect further 
data up to higher transverse momenta and with smaller 
uncertainties in the near future for all the observables 
discussed here. 
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FIG. 4: (Color Online) Pseudorapidity asymmetry, Yaaym for different hadrons at \ri\ < 0.5. The solid line represents the EKS 
nPDFs while the dashed line is obtained with the FGS nPDFs. The dot-dashed line corresponds to the HKN nPDFs, and filled 
triangles denote the STAR data. The solid error bars are systematic errors, while dashed error bars represent statistical errors. 
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FIG. 5: (Color Online) Pseudorapidity asymmetry, Yaaym for different hadrons at 0.5 < |?7| < 1.0. The solid line uses EKS 
nPDFs, while the dashed line is obtained with the FGS nPDFs. The dot-dashed line corresponds to the HKN nPDFs, and 
filled triangles denote the STAR data. The solid error bars are systematic errors, while dashed error bars represent statistical 
errors. 



